Inverse Spin Hall Effect in SNS Josephson Junctions 
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We consider DC supercurrents in SNS junctions. Spin-orbit coupling in combination with Zeeman 
fields can induce an effective vector potential in the normal conductor. As a consequence, an out- 
of-plane spin-density varying along the transverse direction causes a longitudinal phase difference 
between the superconducting terminals. The resulting equilibrium phase coherent supercurrent is 
analogue to the non-equilibrium inverse spin Hall effect in normal conductors. We explicitly compute 
the effect for the Rashba spin orbit coupling in a disordered two-dimensional electron gas with an 
inhomogeneous perpendicular Zeeman field. 

PACS numbers: 



The spin-Hall effect (SHE) and inverse SHE (ISHE) 
are remarkable demonstrations of the influence of the 
spin-orbit coupling on electron transport. Via this cou- 
pling, a longitudinal electric current can induce a perpen- 
dicular spin current and vice versa. These effects take 
place in metals and semiconductors, where the spin-orbit 
interaction (SOI) arises from impurity scattering [l[ or 
band structure effects Utilizing spin injection, SHE, 
and ISHE, electron spins can be controlled, as recently 
demonstrated experimentally Q. 

We discuss the intrinsic SHE and ISHE, where the 
dominant spin-orbit coupling is from the electron band 
structure. The study of SHE has been focused on nor- 
mal conductors, e.g. normal metals and semiconductors. 
New, interesting, and rich physics occurs in supercon- 
ductors where electron transport is dissipationless and 
the ground state exhibits macroscopic coherence. Some 
superconductivity induced features of the intrinsic SHE 
has recently been analyzed in bulk superconductors Q 
and SNS Josephson junctions Q. The latter work re- 
vealed an equilibrium spin accumulation at lateral sam- 
ple edges, similar to non-equilibrium spin accumulation 
in normal conductors, but the spin Hall current vanished 
due to time-reversal symmetry in the DC Josephson ef- 
fect. 

We focus on ISHE in Josephson junctions. There are 
two scenarios depending on how the spin current (den- 
sity) is created in the normal metal. In a dissipative 
setup, additional normal/ferromagnetic terminals in the 
transverse direction inject a non-equilibrium spin cur- 
rent. Subsequently, the ISHE induces an electric poten- 
tial difference Vsh between superconducting terminals, 
causing Josephson oscillations at frequency 2eVsH/fi- 
Transport is dissipative due to the spin flow between 
transverse normal/ferromagnetic terminals. This phe- 
nomenon is interesting from an experimental point of 
view and we will study it quantitatively elsewhere, but 
we consider here a dissipationless effect. 

We predict a novel inverse dissipationless SHE: An 
out-of-plane equilibrium spin density spatially varying in 
the transverse direction induces a longitudinal electric 



supercurrent. Equivalently, it induces a phase shift be- 
tween two superconducting terminals. In general, since 
the equilibrium spin-density controls ISHE, Zeeman in- 
teraction from magnetic or exchange fields manipulates 
the resulting Josephson supercurrent. As an explicit il- 
lustration, we consider the interplay of spin-orbit cou- 
pling and Zeeman fields in a disordered two-dimensional 
electron gas (2DEG), and compute the magnitude of the 
equilibrium Josephson ISHE. 

The interplay of Zeeman field and SOI leading to an 
effective phase difference between superconducting termi- 
nals has recently also been studied in two quite different 
systems, but neither exhibits the ISHE we discuss: A 
supercurrent in response to a spatially homogenous mag- 
netic field has been predicted for Josephson tunneling 
through a ID wire Q and appears in numerical simula- 
tions of the superconducting transport through a ballis- 
tic point contact 0] in a spatially homogenous parallel 
magnetic field. In additon to our main finding of an 
inverse SHE, we provide an improved understanding of 
these phenomena by showing how the interplay of Zee- 
man field and SOI can result in the appearance of an 
effective electromagnetic vector potential. Such a vec- 
tor potential, in direct analogy with the Meissner effect, 
gives rise to a supercurrent. 

Let us outline our model. The spin-orbit interaction 
arises from the band structure, Hso = cr • h^, where 
cr = ((T^, (Jy,(Jz) is a vector of Pauli matrices. We ussume 
that the spin-orbit field hfc is given by Rashba SOI where 
hx = aky and hy ~ —akx- Two examples of spin density 
manipulations in 2DEG will be considered: i) a perpen- 
dicular to 2DEG Zeeman field spatially varying in the 
transverse direction y, as shown in Fig. 1 and ii) homo- 
geneous Zeeman field directed along the y-axis. We will 
show that setup i) exhibits an equilibrium inverse SHE. 
Setup ii) also changes the current-phase relation in SNS 
contacts. All relevant length scales are assumed larger 
than the mean free path I = vpT, and we are in the metal- 
lic regime kpl 3> 1, where hp and vp are the Fermi wave- 
vector and velocity, respectively. These conditions allow 
a diffusion approximation in the description of electron 
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FIG. 1: A SNS Josephson junction. Interplay between Rashba 
spin-orbit interaction and Zeeman splitting in a normal 2D 
film (N) induces a phase difference between order parameters 
of two superconducting terminals (S). An inhomogenous Zee- 
man interaction can be created by e.g. a ferromagnetic layer 
on top of the film or magnetic impurities. Another possible 
configuration (not shown) is a uniform field parallel to the 
y-axis 



transport. In this regime, the transport properties are de- 
scribed by a generahzed Usadel equation, which we will 
now derive. The resulting Usadel equation is similar to 
the one in Ref. Q , but important non-trivial new terms 
essential for the effects we discuss are added due to the 
Zeeman interaction Hz{v) = azHz{r) + c7yHy{r), where 
Hz {Hy) are the perpendicular (in-plane) components of 
the Zeeman field. We start from the anomalous retarded 
thermal equilibrium Green function i^Q^(r, k, a;), which 
is the Fourier transform of 

Fa/3 = -l{[^a.{v+^,t),i;p{v^^,t'] + )e{t-t') (1) 

with respect to the relative coordinate p and relative time 
t ^ t' . It is convenient to use a singlet-triplet basis rep- 
resenting the Green function. 



F. 



where (3 denotes a spin projection opposite to /3, = 
a^p ± icT^^p. Fs denotes the singlet component. Fq and 
F±i are triplet components corresponding to and ±1 
projections of the Cooper's pair total spin on the z-axis. 
Using a standard method starting from Gor'kov equa- 
tions [lol | , we derive the semiclassical equation 



2ttN, 



-[G?ivi/ + *G°2 



(3) 



where subscripts n and m attain the values 0,±1, or s, 
and 



K = 2tj - vq - 2Jhk - S - B. 



Here q = — «V, J is the 3x3 matrix spin 1 operator in the 
triplet subspace, and operators S and B provide mixing 
of triplet and singlet components: 



q dh^ 

■ Bqs — Bsn — 2Hz 



V2 dk 

iV2Hy, 



(5) 



where = h'^ i: ih^ . In the right-hand side of Eq. ([3]) 
^' = F and the unperturbed retarded Green func- 
tions are 
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{ujTEk~(T-W- <7,Hz =F cTyHy + iT) 



(6) 

The diffusion equation can be derived from Eq. ([3]) by 
expanding the operator (1 — iTK)~^ for small tK and 
averaging over k. The resulting Usadel equation is 

2iuj-^ = t{ (^-iv • ^ + 2J • hfe^ )f«' - A/^' , (7) 

where the angular brackets denote averaging over the 
Fermi surface. The matrix M originates from the SOI 
and the Zeeman interaction expressed via the operators 
S and B. Its off-diagonal terms describe singlet-triplet 
transitions. The relevant matrix elements for our further 
analysis arc : 



V2 



Ajj-'i 1 

Mos = -2iH, , 



(8) 



where ~ q'^dh^/dk'' and = h^{\r ■ q) so that e.g. 
the singlet-singlet diagonal clement is proportional to V^: . 

In order to understand some of the underlying physics 
described by Eq. we will demonstrate that SOI in 
combination with the Zeeman field gives rise to an ef- 
fective Meissner effect. Let us first discuss this in the 
most transparent "local" approximation when the SOI is 
strong enough/the system long enough, so that the spin 
diffusion length Lso = vp j ^ D jV so ^ F, ^Jd/T, where 
F is the length of the junction, Tso = 2r(/i^) p is the spin 
relaxation rate and D = vj^T /2 is the diffusion constant. 
In this approximation derivatives in triplet parts of Eq. 
(O can be disregarded, except in a narrow range ~ Lso 
near the boundaries. is assumed to vary slowly on 
the Fgo scale. Expressing the triplet components of ^' 
via the singlet g and substituting them into the singlet 
projection of Eq.©, the latter takes the form 
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ox 



(9) 



(4) where A is a real coefficient obtained from the equation 



2iA\',: = Mss 



1 



■ E 

m=±l 



MsmM„ 



(10) 



Here we have only included dominant terms proportional 
to o^dHz/dy and aHy, where a is the SOI coupHng con- 
stant. Higher order contributions to ^ proportional to 
and have been disregarded. 
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The diffusion equation ([9]) demonstrates that cAje 
is an effective weak electromagnetic vector-potential. 
Therefore, similar to the Meissner effect it will induce 
a supercurrent. To order the solution of Eq. is 
= ^'^exp(ia;A/D), where 4"° satisfies Eq. © with 
A = 0. The exponential factor gives rise to an additional 
phase difference Q = LA/D between the superconducting 
terminals, the Josephson current is jcSin(0 + 9), where 
(j) is the initial phase difference between the terminals 
and jc is the critical current determined by the function 
The coefficient A is simple for Rashba SOI. For a 
parallel Zeeman field A ~ AarHy. For a perpendicular 
field it vanishes, which is expected since it is similar to 
the behavior of the spin Hall conductance. Continuing 
such an analogy, one can expect that A for the cubic 



Dresselhaus [JJl SOI [ll|. 



In order to find a finite ISHE even for the Rashba SOI, 
we must extend our consideration beyond the local ap- 
proximation. In this case the diffusion equation ([7]) can- 
not be reduced to the simple form ([9]). We consider su- 
perconducting leads with equal real order parameters A 
connected via two SN interfaces wit a low transparency 
t. The barriers are assumed to extend into the 2DEG 
under the superconducting leads, so that the range of a 
free electron motion is between xl and xr at the left and 
right leads, respectively. Depending on contact fabrica- 
tion, other models can be similarly studied. For exam- 
ple, the electrons in the 2DEG could move freely under 
contacts, with the barriers present only in z-direction, as 
shown in Fig.l. The choice of the model is not important 
for the main qualitative results obtained below. 

To the lowest order in the tunneling transparency t, 
the superconducting current can be expressed as a 
sum over Matsubara frequencies lo = 7r(2n + 1)T: 



ireT 
2RINf 



A2 



Im[/ rfydy7,,(ri,4)], (11) 



Eq. ^ with hard- wall boundary conditions, Vxfsf' 
'At X — XL and x = xr, while the triplet components in 
the case of Rashba SOI satisfy the boundary condition 
{iL,oV., + 2Jy)f ^0 [13. 

To illustrate the ISHE, we consider the case of Rashba 
SOI with finite Hz, Hy = 0, and A ^ w. The parameter 
of interest is the effective phase difference 



7lm[/(iydy7,s(rL,r'fl)] 



A2 



(13) 



Fro m Eq. 0, 6'eff = C0, where C = 2l\'yHzT/kFL 
and VyHz denotes the average value of the the mag- 
netic field gradient in the contact range, (j) is shown in 
Fig. 2 as a function of ratio of the spin relaxation rate 
Tso — 2Ta^kp versus the Thouless energy Et = D/L^, 
a convenient measure of the SOI strength. For large 
SOI the "local" approximation is obtained by using in 
the second term of Eq. the approximate form of 

f±i±i = 2/oo = —S{ri — r2)/Tso- In this case both terms 
in (|12p are proportional to and precisely cancel each 
other, as in the factor A in Beyond this leading 

"local" approximation there are terms increasing slower 
than a^. They contribute to Fig. 2. 

Larger Zeeman fields cannot be treated perturbatively. 
A strong depairing effect takes place when the character- 
istic length Lz = y^D/2H is small, Lz ^ min(i, L^o). 
Then, for H — H^, both 'fs and 'I'o decay exponen- 
tially near contacts with superconducting terminals and 
the latter become effectively disconnected. On the other 
hand, as it follows from Eq. '^±i components are 

not subject to the depairing effect and can propagate at 
the relatively large distance ~ Lso- Such a long-range 
triplet effect has been studied in SFS junctions, where 
a link between triplet and singlet Cooper pairs has been 
induced by an inhomogeneous (rotating) magnetization 



where Ri, is the boundary resistance [14| |. v^^r 
{xL/R,y) and /ab(ri,r'^), with a,b = 0, ±l,s, is the 
Green function of Eq. ([7]), i.e. a solution of Eq. ([7|) 
with a delta source in its right-hand side. The equations 
for retarded and advanced functions must be properly 
continued to the upper and lower complex semiplanes of 
uj, respectively. Treating M in ([7|) perturbatively one can 
express the correction to /is°''(ri, r^) as 

SUrL^r'n) = -J dr/i«) (r^, r)Af,,/(o)(r, r^) + 
J dridr2/i°)(rL,ri)Af,„ x 

mm' 

/i°i'(ri,r2)M„,Ji°)(r2,ri,), (12) 

where the unperturbed diffusion propagators /is°''(r, r') 
and /^^/ (r, r')i with m, m' = or ±1, are obtained from 








FIG. 2: The phase difference versus a ratio of the spin re- 
laxation rate and the Thouless energy, at 0.5< ksT / Et < 8. 
The parameter C is described in the text. 
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(see [lO[ and references therein). In our case, a coupling 
of ^±1 to \E'o and ^'s can be provided by SOI through 
the matrix elements M±is and the spin precession opera- 
tor i?±i,o = -«4r((J±i,o • hk)iv ■ -§p))f originating from 
the first term in the right-hand side of Eq. ([7]). Indeed, 
assuming that ^ Fgo and Et, it is easy to show that 
the modified Eq. p2p is represented by its second term, 
where the integrand has the form 

/io"^ {xL,x)Ro.nfiSL{x, x')M^J^\x\ xr). (14) 

The unperturbed functions /^q^ and /is"'' arc obtained 
from s and projections of Eq. where the preces- 
sion term and all My , except Mso and Mqs are ignored. 
The physics of the process described by (|14p is clear: the 
magnetic field mixes 0-triplet and singlet components of 
the pairing function within the short range near the left 
boundary. Further, due to the spin precession in the SOI 
field the 0-triplet transforms to ±1 triplet components. 
The latter propagate to the right contact where they con- 
vert to the singlet through M±is. Integrating ([T^ over x 
and x' gives a power law dependence of \m\8fss{xL, xr)] 
on the magnetic field: 

Im[5UxL,XR)] (X {\H,{y^)\--'/^-\H,{y^)\-^/^), (15) 

where y\ and j/2 are y-coordinates of the junction edges. 

In contrast to a perpendicular Zeeman field, in a paral- 
lel field ±1 triplets exponentially decay near boundaries, 
as can be seen from Eqs. (|7l8p . So they cannot provide 
a long-range link between superconducting terminals. 

In conclusion, an analogue to the ISHE exists in DC 
Josephson SNS junctions. Unlike the normal ISHE, the 
supercurrent through the SNS contact can be induced 
by a static Zeeman interaction by magnetic or exchange 
fields oriented normal to the 2DEG and varying in the 
direction transverse to the electric current. A destructive 
depairing effect of the strong Zeeman field is diminished 
by Rashba SOI leading to a power-low dependence on this 
field. We show that a supercurrent through the junction 
can also be induced by a uniform parallel Zeeman field, 
corroborating thus the numerical analysis of Ref. 
On the other hand, the depairing effect of such a field 
was found to be strong (exponential). In both cases an 
appearance of the supercurrent can be explained in terms 
of the Meissner effect produced by an effective vector 
potential, which is a combined effect of the Zeeman field 
and Rashba spin-orbit interaction. 

We considered the diffusive transport regime which is 
relevant in low mobility metals and (magnetic) semicon- 



ductors. Furthermore, the diffusive regime, allows an 
elucidation of the main physics and parameters govern- 
ing this phenomena. We expect a strong Josephson ISHE 
in ballistic junctions containing a metallic normal layer 
with a strong Rashba interaction, for example in Bi films 
on some substrates [l^ . Ballistic quantum wells of nar- 
row gap semiconductor are also expected to exhibit an 
increased Josephson ISHE. 
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